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INTRODUCTIC»l 


This research program, covering the 53-month period, February 1978 
through June 1981, and performed under NASA Grant NSG-7418, has been a 
continuation of the original il^llo lunar passive seismic experiment 
(NASA Contracts NAS 9-5957, NAS 9-13149, NAS 9-14581 and NAS 9-14810) 
following the termination of the data acquisition phase of the experiment. 
The principal objectives of the program have been (1) to complete process- 
ing of the acquired data to provide lunar seismic data flies accessible 
to all potential Investigators, and (2) to analyze the data for timely 
dissemination of significant new results to the lunar and planet 2 ury 
science community. 

Thlis final report Is a suitanary of acconpllshroent during the period 
of the Grant. Detailed results of the Investigation have been published 
during the course of the study. These publications are listed In the 
PUBLICATIONS section of this report. Abstracts of these publications are 
includ-id In the Appendix. 


SUMMARY OF RESULTS 

Task 1. Data processing, event identification and cataloguing 

The scientific data transmitted continuously from all ALSEP (Apollo 
Lunar Surface Experiment Package) stations on the moon and recorded on 
instrumentation tapes at receiving stations distributed around the earth 
were sent to our Laboratory in Galveston for processing. Tho processing 
produced sets ot con^uter-con^iatible digital tapes, from which various 
other data sets convenient for analysis were generated. We read the 
seismograms produced, identified and classified various types of seismic 
events, emd catalogued the detected events. In all we catalogued 12,558 
events detected on the long-period seismograms. Much more numerous events 
detected only on short period selsnx>grams, representing seismic events of 
local nature, were not catalogued. The details of the data processing are 
given in Nakamura et [19801*. The processed data now availedjle aure 
listed below. The following codes are used In parentheses to Indicate 
the location of each data set: 

N: National Space Science Data Center, NASA Goddard Space 

Flight Center, Code 601, Greenbelt, MD 20771. 

G; Galveston Marine Geophysics Laboratory, Institute for 
Geophysics, The University of Texas at Austin, 

Galveston, TX 77550-2768 


•All references are listed in PUBLICATIONS section. 
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1. Continuous AI>SEP normal-blt-rate data on magnetic tapes (N,G) . 

This set of 1151 reels of 9-track magnetic tapes and 3490 reels of 
7-track magnetic tapes (6 only) contains continuous digital data for 
the time period 1 March 1976 through 30 September 1977 for all stations. 
The data include all PSE (Passive Seismic Eiqperiment) seismic data 
as well as those of other experiments except LSPE (Lunar Seismic 
Profiling Experiment) . (The data before 1 March 1976 were processed 
at the Johnson Space Center , and are available from NSSDC.) 

2. Continuous ALSEP high-bit-rate data on magnetic tapes (N,G) . 

This set of 336 reels of 9-track magnetic tapes contains continuous 
digital data of LSPE in the listening mode for the time period 15 
August 1976 through 25 April 1977. 

3. Continuous, compressed-time-scale (75 rm/hr) PSE seismograms on micro- 
films (M,G). 

This data set continuously covers the entire dxaration of the 
experiment from 21 July 1969 to 30 September 1977. (Originals on 
white paper cure in Galveston.) 

4. Continuous, compressed-time-scale (90 mm/hr) LSPE seismograms on micor- 
films (M,G) . 

This data set continuously covers the duration from 15 August 1976 
to 25 April 1977. (Originals on white paper are in Galveston.) 

5. Long-period event catalogs (N,G) . 

These catalogs list all seismic events visually detected on the 
long-period compressed-time-scale seismograms. The catalogs con- 
tain signal start and stop times to the nearest minute, an^litudes 
as measured on seismograms, availability of expanded-time-scale 
seismograms, quality of data, event classification, and comments. 

They were initially pxiblished in nine volumes, first by Apollo 
missions amd later by the year, as listed in PUBLICATIONS section. 

The final version contains all 12,558 detected events in a single 
volume. They are available on print, microfiches, or a magnetic 
tape. 

6. Event tapes (N,G) . 

This set of 1172 reels of 7-track magnetic tapes contains digital 
data for all long-period seismic events listed in the event cata- 
log except those detected only at Apollo 16 station with measured 
amplitude less than 5 mm. 

7. Compressed-tiine-scale (75 mm/hr) seismograms of events on microfilms 
(N,G) . 

This data set contains compressed seismograims of all events recorded 
on the event tapes. (Originals on white paper are in Galveston.) 

8. Expanded-time-scale (100 mm/min) seismogrcuns of events (G) . 

This data set, produced on translucent papers, contains expanded 
seismograms for the first 10 minutes of all events recorded on 
the event tapes. 
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9. Artifical Impact data on magnetic tapes (N,G) 

This set of 5 reels of 7-track tapes contain digital seismic data 
for all 9 artificial inpacts. A 9-track version is also available 
on a single reel (G only) . 

10. Major meteoroid inpact data on magnetic tapes (N,G) 

This set of 32 reels of 7-track tapes contains digital seismic data 
for 98 meteoroid Inpact events which registered 10 mm or more on 
our standard conpressed seismograms at two or more stations. A 
9-track version is also available as a set of 8 reels of tapes (G 
only) . 

11. Shallow moonqueJce data on magnetic tapes (N,G) 

This set of 9 reelr of 7-trcck tapes contains digital seismic data 
for all 28 rha" Lo- moonq[ueJces detected. A 9-track version is also 
available as a set of 2 reels of tapes (G only). 

12. Major deep moonquake data on magnetic tapes (G) 

This set of 196 reels of 7-track tapes contains digital seismic 
data for initial 45 minute sections of signals from 1074 deep moon- 
quakes at 33 source regions from which 10 or more moonquakes were 
detected. A highly condensed version containing only the :.nitial 
8 minute sections of long-period signals is also available on a 
single reel of 9-track tape. 

Not listed cibove, but available in Galveston are a set of 7267 reels 
of 7-track magnetic tapes that contain continuous PSE digital data for the 
period 21 July 1969 through 29 February 1976, produced at the Johnson Space 
Center, and various data sets produced for our own analysis use. The latter 
include a set of stacked seismogreuns for 33 major deep moonquake source 
regions, stacked digital seismic data on a 5-track tape, a set of 24 reels 
of 7-track magnetic tapes containing Fourier-trcinsformed seismic data for 
major impacts and moonquakes, and 30-second low-pass filtered continuous 
long-period seismic data for the time period 1 March 1976 through 27 March 
1977. 


Task 2. Data analysis 
Seismicity 


We made a detailed amalysis of the bovirce mechamism of deep moonquakes 
at Aj source region. Since various factors, including paucity of seismic 
stations and complexity of lunar seismic signals, precluded use of conven- 
tional techniques such as distribution of initial motions, we used cross- 
phase spectra and amplitude distributions [Nakaunura, 1978] amd polarization 
of aurrivals [Koyaria and Nakamura, 1980] to study the source mechanisms. 

The results show that all mconqiiake foci within the A^ source region are 
distributed in a nearly horizontal plauiar regions less than 1 km in diameter, 
and that shear fractures occur along a horizontal plane with the slip 


' k 
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direction rotating with the shifting tidal stress field. Thus we conclude 
that these deep moonquakes represent a process of single storage and release 
of tidal energy rather than a release of accumulated tectonic stress. 

We have identified 107 distinct source regions of deep moonquakes > and 
located 52 of them using our latest seismic velocity model of the lunar 
interior [Nakamura ^ al. 1982] . Several linear patterns of eplcentral 
distribution are recognized. Some of them appear to reflect the bilateral 
symmetry of the tidal force which produce these moonquedces. There appezu: 
to be two zones where many deep moonquakes are concentrated: one is just 

below 800 km depth, and the other is just above 1000 km depth. These two 
levels may suggest internal boundaries of differing physical properties. 

Shallow moonquakes appear to be the only tectonic quakes in the moon. 

The vaoriation of observed seismic amplitude with distance indicates that 
most, if not all, of them originate in the top part of the liinau: upper mantle 
[Nakamura et 1979] , thus suggesting high accumulation of tectonic stress 
in this region of the lunar interior at present. Shallow moonquadces show 
a remarkable slmilarlly to intraplate earthquakes [Nakamura, 1980]. Al- 
though they do not produce laorge chemges in sxirface geologic featxures like 
earthq\uJces along plate boundaries, these Infrequent shallow moonquakes, 
like intraplate earthquakes, could be very large because a large amoxint of 
strain energy can be accumulated within a lithospheric plate before it is 
released by a quake. The very thick lithosphere of the moon does not allow 
itself to break up, thus leaving only this type of tectonic quakes to occur 
there. 

Meteoroid in^acts detected by the long-period Reismometers of the 
Apollo lunar seismic network reveal the distribution of small interplaneteury 
bodies near the earth-moon orbit in the mass range from about 100 g to 100 
kg. The observed inpacts cue not random, but show clustering. Some, but 
not all, of the clusters appear to be correlated with known meteor streams 
[Dorman ^ 1978] . Further studies of these impacts and much more 

numerous smaller impacts detected by the short-period seismometers are needed 
to clarify the nature of these interplanetary bodies. 


Internal structure 


Seismic data can be used in several different ways to infer the structure 
of planetary interior. Though perhaps the most common way is to use travel 
times of body wave arrivals, such things as dispersion and polarization of 
surface waves and anplitudes of body waves can be used. We noticed that 
observed anplitudes of seismic coda were significantly different for hori- 
zontal and vertical conponents and that they were also frequency dependent. 
Interpreting this to be due to the ellipticity of Rayleigh wave particle 
motion, we derived the ne^Lr-surface velocity structure at the Apollo 12, 

14, 15 and 16 landing sites [Horvath et al . 1980] . At these sites, the 
shear-wave velocity-depth functions are generally represented by a rapid 
increase in the top 10 to 25 m layer followed by a more gradual increase 
to a depth of 150 to 200 m. Details of velocity variations appear to be 
correlated with the history of the lamding sites. 
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Structure of the deep Interior of the moon can be estimated from the 
arrival times of distant seismic events. An inversion of our latest set 
of arrival-time readings show that both P and S waves show a negative ve- 
locity gradient in the upper mantle and an increase of velocity in the 
ndddle mantle between the depths of 500 and 1000 km [Nakamura, 1982] . Such 
velocity variations are consistent with an upper maintle of nearly uniform 
con^sition consisting of a mixture of olivines and pyroxenes, and a middle 
mantle of higher concentration in Mg-rich olivines. The olivine enrich- 
ment may suggest melting of the moon at least to a depth of 1000 km in its 
early history. 

We have attenpted to identify secondary seismic arrivals reflected from 
possible discontinuities in the deep lunar interior using various signal 
enhancement techniques [Horvath, 1981]. The results, however, were incon- 
clusive. Attempts tc detect coherent surface waves and free oscillations 
of the moon have bean unsuccessful. 

Seismic Q in the moon except its deep interior is extremely high. A 
recent study shows that S-wave Q in the lunar upper mantle is frequency de- 
pendent ^md is significantly higher than P-wave Q at high frequencies 
[Nakamura emd Koyama, 1982] . In absence of other dissipation mechanisms 
normally operating in the earth, congressional heat loss may be a dominant 
dissipation mechanism in the Ixinar upper mantle. 

Lunar seismogrcims are very con^lex because of intense scattering of 
seismic waves near the lunar surface. Certain sin^lified mathematical 
models can explain some of the observed characteristics of lunar seismograms 
[Malin and Nakamiira, 1901] , but further studies are needed for a better 
understanding of this phenomenon. 

In summary, the analysis of the lunar seismic data has revealed many 
interesting properties of the Ivinar interior. The moon is clearly 
differentiated, having a distinct crust amd a mantle, and possibly a core. 

The interior of the moon is still of high teit^erature. However, because 
of its great thickness, the lithosphere of the moon is not broken up into 
separate plates like those on the earth. Rare but large tectonic moonquakes, 
similar to intraplate earthquakes, occur in the shallow part of the litho- 
sphere, while numerous but small tidally induced moonquakes occur near the 
bottom of the lithosphere. Overall reviews of the results of the lunar 
seismic experiment are found in Nakeunura [1981] and Nakaimura ejt al. [1982] . 
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Ai moonquakes: Source distribution and mechanism 

Yosio Nakamura 

University of Texas, Marine Science Institute. Geophysics Laboratory, Galveston. Texas 77550 


Abstract — Two novel analysis techniques have been used to study focal mechanisms of Ai moon- 
quakes. which occur at the most active of the 80 deep-moonquakc source regions identified to date. The 
first is the cross-spectral analysis of matching signal waveforms. Cross phase spectra clearly show rela- 
tive polarity of signals and give precise time relationships between signals, thus enabling accurate deter- 
mination of relative source locations within the source region. The second is the analysis of amplitude 
ratios, which reveal the patterns of P- and S-wave radiations around the source region. Ai moonquake 
foci are found to be concentrated in a nearly horizontal planar region less than 1 km in diameter. Their 
focal mechanism may be described by slippages along a horizontal plane with the slip direction rotating 
with and predominantly controlled by the shifting tidal stress field. Thus, Ai moonquakes represent a 
process of simple storage and release of tidal energy with no evidence for a release of accumulated 
tectonic stress. 
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Passive seismic experiment: a summary of current status 

G. V. Latham, H. J. Dorman, P. Horvath, A. K. Ibrahim, J. Koyama, 

and Y. Nakamura 

University of Texas, Marine Science Institute, Geophysics Laboratory, Galveston. Texas 77550 


Abstract — The data set obtained from the 4-station Apollo seismic network is now complete. It 
includes signals from approximately 11,800 events of various types. Four data sets for use by other 
investigators, through the NSSDC, are in preparation. Some refinement of the lunar model based on 
seismic data can be expected, but its gross features remain as presented two years ago. The existence of 
a small, molten core remains dependent upon the analysis of signals from a single, far-side impact. 
Analysis of secondary arrivals from other sources may eventually resolve this issue, as well as continued 
refinement of the magnetic field measurements. Evidence of considerable lateral heterogeneity within 
the moon continues to build. The mystery of the much lower meteoroid flux estimate derived from 
lunar seismic measurements, as compared with earth-based estimates, remains; although, significant 
correlations oetween terrestrial and lunar observations are beginning to emerge. 
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On the time-varying properties of the lunar 
seismic meteoroid population 

J. Dorman, S. Evans, Y. Nakamura, and G. Latham 

University of Texas. Marine Science Institute. Geophysics Laboratory. Galveston, Tsxas 77550 

Abstract— Strong short-term fluctuations of meteoroid impact rates are evident in a list of 1557 events 
derived from Apollo lunar seismic data. Times of fall and seismic signal amplitudes are considered in 
relation to the possible orbits and identification of the impacting objects. The entire lunar surface is the 
effective collector; and the Apollo network data, gathered between 1970 and 1977, reflect the recog- 
nized abundance of large meteorite falls from early April through July. There is also some indication 
that briefly increased counts represent fragments of about 100 g or greater belonging to several meteor 
streams: Quadrantids, Aquarids. Perseids. Oriunids. Leonids. Geminids, and possibly others as welt. 
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MARTIAN WIND ACTIVITY DETECTED 3Y A SEISMOMETER AT VIKING LANDER 2 SITE 

Yosio Nakamura 

Geophysics Laboratory, Marine Science Institute, The University of Texas, Galveston, Texas 77550 

Don L. Anderson 

Seismological Laboratory, California Institute of Technology, Pasadena, California 91125 


Abstract . The seismic 'background noise' 
detected by the seismometer on the Viking lander 
2 has an extremely high correlation with the 
measured wind speed at the landing site. When 
displayed in a compressed form, the nearly 
continuous seismic data clearly exhibit the 
diurnal as well as seasonal variations of the 
Martian wind activity. A preliminary spectral 
analysis of th,' long-term variation of the back- 
ground noise indicates persistent spectral peaks 
at periods near 1.5, 3, 7 and 10-20 sols for the 
first 560 sols (0.84 Martian year) of observa- 
tion. 



16 

CKiGINAL tS 

OF POOR QUALITY 


Fnn t.imur P/i/wl Set lOth (1979). p. :299-:i09 

Printed in ihe United S(ate> ol America 


Shallow moonquakes: Depth, distribution and 
implications as to the present state of the lunar 

interior 

Yosio Nakamura, Gary V. Latham, H. James Dorman, 
Abou-Bakr K. Ibrahim, Junji Koyama, and Peter Horvath 

The University of Texas. Marine Science Institute. Galveston Geophysics Laboratory. 

Galveston. Texas 77550 


Abstract — The observed seismic amplitudes of HFT (high-frequency teleseismic) events do not vary 
with distance as expected for surface sources, but are consistent with sources in the upper mantle of 
the moon. Thus, the upper mantle of the moon is the only zone where tectonic stresses deriving from 
differential thermal contraction and expansion of the lunar interior are presently high enough to cause 
moonquakes. The distribution of shallow moonquake epicenters suggests a possible correlation with 
impact basins, implying a lasting tectonic influence of impact basins long after their formation. The 
finite depths now assigned to these shallow moonquakes necessitate further revision to the seismic 
structural model of the lunar interior. 
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HOW WE PROCESSED APOLLO LUNAR SEISMIC DATA • 

YOSIO NAKAMURA, GARY V. LATHAM and H. JAMES DORMAN 

University of Texas Marine Science Institute, Galveston Geophysics Laboratory, Galveston, TX 77550 (U.S.A.) 

The Apollo lunar seismic station network gathered data continuously at a rate of 3 X 10* bits pet day for nearly 
eight years until the termination in September, 1977. The data were processed and analyzed using a PDP-15 mini- 
computer. On average. 1500 long-period seismic events were detected yearly. Automatic event detection and identi- 
fication schemes proved unsuccessful because of occasional high noise levels and, above all, the risk of overlooking 
unusual natural events. The processing procedures which were finally chosen consist of plotting all the data on a 
compressed lime scale, visually picking events from the plots, transferring event data to separate sets of tapes and 
performing detailed analyses using the latter. Many problems remain, especially in the automatic processing of extra- 
terrestrial seismic signals. 
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Lunar Near-Surface Shear Wave Velocities at the Apollo Landing Sites 
us Inferred From Spectral Amplitude Ratios 

Peter Horvath,' Gary V. Latham, Yosio Nakamura, and H. James Dorman 

University of Texas, Marine Science Institute, Gelveslon Geophysics Laboratory, Galveston, Texas 77550 

We reexamined the horizonul-to-vertical amplitude ratios of the long-period seismograms to deter- 
mine the shear wave velocity distributions at the Apollo 12, 14. 15, and 16 lunar landing sites. Average 
spectral ratios, computed from a number of impact signals, were compared with spectral ratios calculated 
for the fundamental mode Rayleigh waves in media consisting of homogereous, isotropic, horizontal lay- 
ers. The shear velocities of the best fitting models at the different sites resemble each other and differ 
from the average for all sites by not more than 20% except for the bottom layer at station 14. The shear 
velocities increase from 40 m/s at the surface to about 400 m/s at depths between 95 and 160 m at the 
various sites. Within this depth range the velocity-depth functions are well represented by two piecewise 
linear segments, although the presence of first-order discontinuities cannot be ruled out. 
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Shallow moonquakes: How they compare with earthquakes 


Yosio Nakamura 

The University of Texas. Marine Science Institute, Galveston Geophysics Laboratory. 

Galveston. Texas 77550 


Abstract — Of three types of moonquakes strong enough to be detectable at large distances — deep 
moonquakes, meteoroid impacts and shallow miHinquakes — only shallow moonquakes are similar in 
nature to earthquake.s. A comparison of various characteristics of moonquakes with those ot earth- 
quakes indeed shows a remarkable similarity between shallow moonquakes and iniraplate earth- 
quakes: ( I) their occurrences are not controlled by tides; t2) they appear to occur in locations where 
there is evidence of structural weaknesses; (.1) the relative abundances of small and large quakes (b- 
valucs) are similar, suggesting similar mechanisms; and (4) even the levels of activity may be clo.se. 
The shallow moom,uakes may be quite comparable in nature to intraplate earthquakes, and they may 
be of similai origin. 
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Focal mechanism of deep moonquakes 

Junji Koyama* and Yosio Nakamura 

Marine Science Institute, Geophysics Laboratory, University ofTexas, Galveston. Texas 77550 


Abstract— To elucidate the focal mechanism of deep moonquakes, we analyzed S-wave polarizations 
of deep moonquake signals from the A, source region. At station 12, where the available data are of 
the highest quality, the variation of polarization angle with the anomalistic phase of the moon indicates 
agreement with that expected from the focal mechanism model of Nakamura. Focal mechanism 
solutions were derived for eight A| moonquakes assuming that moonquakes are. like earthquakes, 
caused by a shear fracture on a fault plane. The mechanism solutions generally indicate a nearly 
horizontal or almost vertical faulting. The slip directions estimated from the solutions are different 
from one another, again suggesting variation of focal mechanisms as a function of the tidal phase of 
the moon. 
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Correction of lunar seismograms for instrumental and 
near-surface effects and constraints on the 
velocity structure of the lunar interior 

Peter Horvath’^ 

The University of Texas. Marine Science Institute. Galveston Geophysics Laboratory, 

Galveston. Texas 77550 


Abstract — Predictable digitizing errors were removed from the long-period lunar seismograms, and 
inverse filters were applied to correct for the seismometer response, for the coupling of the 
seismometer to the ground, and for the near-surface structural effects. These signal processing steps 
improved correlations among the long-period seismograms from the various recording sites and 
permitted reading the shear wave arrivals with smaller uncertainty than is possible from the original 
seismograms. To explore the constraints that direct compressional (P) and shear (S) arrival times 
could impose on the velocity structure of the lunar interior, the largest natural imp.acts and shallow 
moonquakes were located using two velocity models. The first model consisted of a thicker crust 
(55 km) and higher compressional and shear velocities (8.1 and 4.6 km/sec. respectively) in the upper 
mantle, and the second model had a thinner crust (45 km) and lower upper-mantle velocities (7.7 and 
4.4 km/sec). The P and S arrival times observed fur impacts and shallow moonquakes could not be 
used to distinguish between these two models because the travel time residuals for both models are 
similar and are well within the uncertainty of the observed travel times. The deconvolved seis- 
mograms have also been used to search for secondary arrivals and thereby to obtain additional 
constraints on the velocity structure. The identification of possible secondary arrivals suggests that a 
velocity discontinuity exists in the crust at a depth between 20 and 50 km. ,ind that (he crust is 
thinner than 50 km. Although evidence for a discontinuity in the mantle was inconclusive, dis- 
continuity at a depth of about 200 to 220 km can be inferred from secondary arrivals. 
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Seismic Q of the Lunar Upper Mantle 

Yosio Nakamura and Junji Koyama* 

Instiiutefor Geophysics, The University of Texas at Austin. Galveston, Texas 77530 

We have determined the Trequency dependence of anclastic attenuation for both P and S waves in the 
upper mantle of the moon in the frequency range of roughly 3 to 8 Hz. The method used is a combination 
of a single-station method, in which is obtained as an integral with an unknown constant, and a 
multiple-station method, the result of which is used to determine the integration constant. The Q 
determined for P waves appears to decrease with increasing frequency, though it is not significantly 
diflerent from a constant Q, while the Q for S waves (Q,) increases with frequency, becoming significantly 
greater than Q, at high frequencies. All Q values remain greater than 4000 within this frequency range. 
Q, above 5 Hz is approximately proportional to the 0.7ih power of frequency. The high Q for shear waves 
at high frequencies suggests that compressional heat loss may be a dominant dissipation mechanism. Some 
uncertainties in the absolute Q values remain because the detailed velocity structure of the lunar upper 
mantle is largely unknown. 


6540 Interior of moon 

APOLLO LUNAR SEISMIC EXPERIMENT - FINAL SUMMARY 
Yosio Nakamura (Institute for Geophysics, The 
University of Texas at Austin, Galveston, Texas, 77550- 
2768), Gary V. Latham, and H. James Dorman 

Processing and initial analysis of the entire set of 
Apollo lunar seismic data collected continuously from 
1969 through 1977 have now been completed. Recent 
results include; 1) better defined deep moonquake 
locations, which appear to be bounded rather sharply 
between about 800 km and 1000 km depths with concentra~ 
tions near both boundaries; and 2) middle mantle 
(<v500 to 1000 km depth) seismic velocities of 
Vp ■ 8.3 ± 0.4 km/sec and Vs ■ 4.6 ± 0.2 km/sec, which 
are significantly higher than previous estimates and 
represent an increase of velocities from the upper 
mantle as opposed to a decrease in previous estimates. 
(Lunar seismic data, moonquakes, meteoroids, seismic 
velocities) . 
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Seismic Velocity Structure of the Lunar Mantle 
Yosio Nakeunura 

The University of Texas r Institute for Geophysics 
Galveston Marine Geophysics Laboratory 


Abstract 

I have inverted the recently completed set of seismic arrival times 
from the Apollo lunar seismic network to estimate the average seismic 
velocities in three sections of the lunar mantle: two for the upper mantle 

cind one for the middle mantle. The method used is a variation of the 
linearized least-square inversion where the inversion is accomplished in 
steps. The stimated average velocities in the upper mantle decrease from 
Vp “ 7.74 km/s and Vg » 4.49 km/s in two section 2 d>ove 270 km depth to 
V_ “ 7.46 km/s and V « 4.25 km/s in the section between 270 and 500 km 

r 8 

depth, confirming the earlier finding of negative gradients based on seismic 
aunplitude variations. The average velocities in the middle mantle between 
the depths of 500 km and 1000 km of Vp » 8.26 km/s and Vg ■ 4.65 km/s are 
significantly higher than those in the upper mantle, contradicting earlier 
estimates based on more limited data. The higher velocities may suggest 
initial melting of the moon down to at least 1000 km depth. 


